Local Average and Scaling Property of 1/ f a Random Fields 
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We use the techniques developed in Q to study the local average of random fields with spectral 
density l/f a . We study their scaling properties and show that the self-similarity of 1// random 
fields is preserved under the local average. We study how the original random fields can be recovered 
from the locally averaged ones. We also study the derivative of the locally averaged random fields as 
a way to get the spectral density. Finally, we propose the generalization of local average by means 
of an arbitrary response function. 

PACS numbers: 02.50.-r,05.40.-a 



INTRODUCTION 

Random fields are associated with probabilities. If 
we want to study a one-dimensional random field Q 
X(t), we would like to specify not only the proba- 
bility density function P(x, t) for this field but also 
the joint probability density functions P{x\,tx;x2,t 2 ),..., 
P{x\,ti;x 2 ,t 2 , —,x n , t n ) and so on. Here the capital let- 
ter X denotes the random field, and the lower case letter 
x denotes the possible value it can take, and t is the 
time. With those probabilities in hand, we can study 
many statistical properties of the random fields. Some 
are first order effects, like the mean or expectation value 

E[X(t)]= J xP(x,t)dx, some are second order effects, 

— oo 

like the 2-point correlation function and covariance, and 
some are higher orders, like the n-point correlation func- 
tions. The n-th order effects are related to the n-point 
joint probability density functions. 

Most random fields have "memories", namely they 
have correlations of their values taken at different times. 
Some don't have correlations and are thus indeed ran- 
dom. Those with non-zero correlations are not that ran- 
dom. Their covariance can be studied as deterministic 
properties. The covariance B(t\,t 2 ) is defined as the 
correlation of the fluctuation around the mean, and is 
associated with the joint probability functions: 
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(1) 



J x 1 x 2 [P(x 1 ,t 1 ;x 2 ,t 2 ) - P{x 1 ,t 1 )P(x 2l t 2 )]dx 1 dx 2 



We can see that the covariance is non-zero only if the 
2-point joint probability density function is not simply 
the product of two probability functions, namely there 
is some entanglement between the two events (x%,ti) 
and (x 2 ,t 2 ) such that the conditional probability den- 
sity function is not the same with the single probabil- 
ity density function: P(x\,ti \ x 2 ,t 2 ) ^ P{x\,t\), and 
P{x 2 ,t 2 | xi, h) y£ P(x 2 ,t 2 ). 



If the random field is stationary, B(ti,t 2 ) is only a 
function of r = t\ — 1 2 : B(ti,t 2 ) — B(t). The covariance 
function B(t) can be more conveniently studied by us- 
ing the Fourier transform. The Fourier transform of the 
covariance function is the spectral density S(u>). With- 
out loss of generality, we can assume that the means of 
the random fields are zero, since the means are determin- 
istic and can be subtracted for simplicity of discussion. 
The random field X(t) can be decomposed into Fourier 
components F(oj): 



X(t) = Re J F(u)t 



(2) 



and thus the covariance of X(t) becomes the Fourier 
transform of the square of F(u) : 



E[X(t 1 )X{t 2 )]=Re i {F 2 {u))e iu ^- t ^dw (3) 



Since F(u>) is the amplitude of each Fourier mode, the 
square of it should be considered as the power or energy 
spectral density: 



(4) 



So we have the Wiener-Khinchinc relation between B{t) 
and S(u)% 
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B(t) — J S{uj) coscurdtu 



(5) 



A class of random fields has the spectral density of the 
form (a is some positive constant): 



S(w) ~ -1. 



(6) 



They are very common in nature, especially for the case 
a = h the 1/f noise, which is reported to be ubiqui- 
tous |2j. The 1/f noise almost present everywhere in 
every kinds of phenomena, from resistors to traffic 
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0, from brain to number theory 0,0, etc. The micro- 
scopic origin of 1/f noise is widely studied, in e.g. 0,0]- 

In real experiment when we measure the value of a ran- 
dom field, we often measure an average value in a small 
domain. The idea of the local average of the random 
fields is thus very practical and important. In this paper 
we will study the local average of the random field with 
spectral density l/w Q , and their scaling properties under 
local average. We will also discuss the derivative and the 
inverse process/ inverse problem of the local average. 

I will first discuss their scaling properties in the co- 
variance functions, both before and after local average. 
Since the spectral density is proportional to ^ , we have 
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B(r) 



— cos uirdu 
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cos u)tcL{u)t). 



So we see 



B(r) 
B(r) 



(for 
(for 



a > 0), 
a = 0) 



(7) 
(8) 
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We can see that the exponent a determines how much 
correlated the ^ random fields are. The bigger the a, 
the more correlated it is. 

If a > 1, B(t) increases with r, it is positively corre- 
lated. For example, in the Brownian motion, we know 
from the Einstein relation that the covariance of the po- 
sition of the Brownian particle increases linearly with 
time (D is the coefficient of diffusion): E[X{ti)X(t2)\ = 
2D(ti — i2), clearly B(t) ~ r, so for Brownian motion, 
a = 2. 

If a < 1, B{t) decreases with r, it is negatively corre- 
lated. If a = 1, it's the 1/f noise, in which B(t) ~ const. 
There is no a frequency scale or time scale in ideal 1/f 
noise, since the covariance function is constant every- 
where, and after a rescaling of the time:r — > nr, it re- 
mains the same: B{nr) — B(t). 

If a = 0, it's the white noise and its spectral density 
S(lu) ~ -4> is constant and B(t) ~ <5(r). 

Generally, the random fields with spectral den- 
sity are self-similar. The curve x vs. t of a realiza- 
tion of these kind of random fields looks like a fractal, 
e.g. 0|"E3 : For general a, under a rescaling:r — > nr, 
B(nr) — n Q_1 _B(r), namely when we magnify r n times 
and simultaneously multiply the random field by a fac- 
tor rS 1 ^ ^/ 2 , the curve x vs. r looks the same. In 
other words, they are fractals with fractal dimensions 
d= 2 + (l-a)/2 H3,[Iil. 

However, these scaling dependence should be appropri- 
ate in the large r regime. In small t, perturbation shows 
B(t) is a Gaussian with maximum B(0) — a x (t). 



Now let's come to the local average of the random fields 
with spectral density -\, 



LOCAL AVERAGE OF l/f a RANDOM FIELDS 

The local average of the random field X(t) with a 
window of length T is defined as 



t+4 



X T (t)= I J Xfadti 



(11) 



There could be many ways to calculate the relations of 
the covariance functions and spectral densities between 
the averaged and unaveraged random fields. 

One way is by defining a variance function j(T) = 

a 2 ■— I 

• It is shown [l| that j(T) is a response to the co- 
variance function: 



7 (T) ~ 1 1(1 - ^)B{r)dT 



(12) 



we have 



For S(oj) ~ ^7, if we use © B{t) ~ r° 

7(T) - r- 1 ." 

Another way is by analyzing the spectral density. It is 
shown in 1] that: 



S Xt M 



sin^T/2) 2 
1 uT/2 J 



S x {u). 



So we have 

Bx T (r) 



-Hoc 



r si n(a/r/2) 2 1 

LOT/2 1 LU a 



cos lutcIuj. 



(13) 



(14) 



When T is very small it's clear that t ™^^ ] 2 ~ 1, 
so both the spectral density and covariance function be- 
fore and after the local average are almost the same. 
Sx T {u) ~ Sx{oj), Bx t {t) « B(t). It remains a ^ spec- 
tral density when T is small. 

When T is not small, consider the function: 



1 



1 



1 



sin (loT/2) coslut = -[coswt — cosa>(r+T) — coso;(r— T)] 



(15) 



So we have 



Bx T {r) 
1 f + °° 

— - / [2 cos ujt — coster + T) 
T Jo 



■ cos U)(T-T)}- 



rduj, 



rp 2 



[2r 



a+l 



(r + TY 



■T\ 



(16) 
(17) 
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When a = 1,Bx t (t) ~ const. The local average of 
1/lu noise is still a 1/oj noise. This is the manifestation 
of the self-similarity of 1 juj noise. 

For a ^ l,when T is small, we can expand T around r 
in the above function, Bx t {t)^ t q_1 , which is the same 
as the unaveraged one. 

If T is big, we can expand r around T, then 



Bx T (r) ~ ^[2r« +1 



2T 



a+l 



a(a + lyT*- 1 }. (18) 



On another hand , we have S x (ui) — uo 2 Sx T (u), so 
we get 



Sx T M 



sin(wT/2) 2 



uT/2 

which is in agreement with eqn (|13|) . 



Sx(w) 



(24) 



INVERSE PROCESS OF LOCAL AVERAGE 



DERIVATIVES OF LOCALLY AVERAGED 
RANDOM FIELDS 

The spectral density of the local averaged random field 
can be achieved in another way in which we first calculate 
the spectral density of the derivative of the local averaged 
field and then use the relation S- (uj) = uj 2 Sx^iu). 

AT 

The derivative of the local averaged random field X? = 

f+i 

^ J X^dh is: 
t-S 



*x = ^= (19) 



The covariance of Xt can be known from the informa- 
tion of the covariance of X, since 

E[X T (t)X T (t + T)} 
= ^{E[X(t+^)X(t+^+r)} 

+E[X(t-l)X(t-^+r)} 
-E[X(t-~)X(t+~+T)] 
-E[X(t+l)X(t-^+r)}} (20) 
If X(t) is stationary, the above equation amounts to 
B Xt (t) - ^ [2B X (r) - B x (r + T) - B x (r - T)] (21) 

If we use B(t) ~ t"" 1 , for a > 0, then 
B^ (t) ~ ^pr"- 1 - (r + T)"" 1 - (|r - T|) Q_1 ] (22) 

It's interesting to compare the covariance function of 
the local average of the random fields before and after 
taking a derivative, i.e. eqn (|18|) and l|22|) . 

Making a Fourier transform of both sides of , 
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[2Sx(w) 



One may ask the question whether we can get the in- 
formation of the original random field X(t) from the in- 
formation of the local averaged one Xx(t). If we have 
two locally averaged fields with very close lengths of the 
window, ie, we measure XT ± (t) and XT 2 (t) in which T\ 
and T 2 are very close to each other (suppose T\ > T 2 ), 
since Xt ± {t) and Xt 2 (t) are integral over range T\ and T2 , 
their difference reveals the information of X^ Tl _ T2 ^(t). 



T 1 X Tl {t)-T 2 X T ^t) 



Ti-T 2 



X 



(Ti-T 2 ) 



Ti+T 2x 



+X {Tl ^ T2) (t+^±^)} 



(25) 



By making the auto-covariance of both sides, we can 
get the cross-covariance of Xt x (t) and Xt 2 (t) . 

If T\ and T 2 are so close that we can make a derivative 
with respect to T, the above equation becomes: 

XT{t)+T^P- = \[X{t-^)+X{t+^)]. (26) 

Together with eqn (|19f) . we may solve X(t) in terms of 
X T (t) : 

x(f) _ x u , T ggrft+f) T dX T jt+%) 
X(t)~X T (t + -)+T — 

(27) 



X(t)-X T (t--)+T — + - 

(28) 

This recovers the original random field X(t) from 
the locally-averaged one Xt(£). The two expressions 
()27J) . (|28|I are consistent. The derivative with respect to 
T is understood as to compare two local averages whose 
lengths of windows are very close to each other. 

The derivative of a local average is the local average of 
the derivative: 



S x (uj) sin 2 (ujT/ 2) (23) 



1 r t+ ~ 

Tj t T 



X{t 1 )dt l = (x(t)) T (30) 
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GENERALIZED LOCAL AVERAGE AND 
RESPONSE FUNCTION 

Local average defined in 1] can be generalized to 
weighted local average, characterized by a normalized 
function h(t): 



X h (t) 



X(tx)h{t x - t)dti 



(31) 



where /_ °° h(t\ — t)dt\ = 1. 

In this language, the ordinary local average CH in 
is the case where the function h(ti — t) is a step function: 



*) 
*) 



1 

T ' 
0, 



for 
for 



1*1 



t € 



T 

2' 



l*i 



t > 



T 



(32) 
(33) 



The generalization of h(t) to other functions might be 
useful when the actual random field we measure is an 
uneven average, namely it is possible that in the averag- 
ing process when we measure in the real experiments, the 
values near the center has bigger weights than those away 
from the center, or inversely. The general h respects this 
experimental detail. For example, perhaps h{t\—t) could 
be a Gaussian centered at t. 

In this case, the spectral density of the weighted local 
average is 



(34) 



where H(u) is the Fourier transform of the function 
h(tx - t). 

We see that this actually means that the local average 
is a response to the original unaveraged random field by 
the response function h. 

I am very grateful to Prof. VanMarcke for helpful dis- 
cussions. 
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